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Abstract 

We present the one-loop QCD corrections to the heUcity amplitudes for the processes 
qq W^W~, Z Z, W^Z, W^j, or Z7, including the subsequent decay of each 
massive vector boson into a pair of leptons. We also give the corresponding tree- 
level amplitudes with an additional gluon radiated off the quark line. Together, these 
amplitudes provide all the necessary input for the calculation of the next-to-leading 
order QCD corrections to the production of any electroweak vector boson pair at hadron 
colliders, including the full spin and decay angle correlations. 
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1 Introduction 



With the increasing energy available at hadron colliders, electroweak gauge boson pair pro- 
duction becomes more and more important. Both experimental collaborations at the Teva- 
tron, CDF [0 and DO p[, have performed studies of pair production processes such as 
pp —>■ W^W~ , ZZ,W^Z,W^'~^, or Z7. These studies considered purely leptonic decays of 
the massive vector bosons in the pair, as well as decays into two jets plus leptons. Already, 
some events have been found above the background, in accordance with Standard Model 
predictions. The amount of available data will increase by roughly a factor of 20 at the up- 
graded Tevatron, and by a factor of 1000 once the Large Hadron Collider at CERN (LHC) 
starts operating. A summary of the present experimental situation and a more complete list 
of earlier studies can be found in 0. 

These processes are quite interesting in several respects. Most of all, they can be used to 
measure the vector boson trilinear couplings predicted by the Standard Model Any kind 
of anomalous couplings, or decays of new particles into vector boson pairs, would result in 
deviations from these predictions. In particular, if the Higgs boson is heavy enough it will 
decay mainly into W^W~ and ZZ pairs Thus, the search for Higgs bosons in the few 
hundred GeV mass range is intimately connected to pair production of vector bosons. A 
detailed understanding of these Standard Model processes is therefore mandatory. 

Due to its importance, hadronic pair production of electroweak vector bosons has received 
a lot of attention in the literature. The tree-level cross-sections for the hadronic production 
of W~^W~, ZZ, W^Z, as well as W'^'j and Zj pairs were computed long ago 0. The 
one-loop {0{as)) QCD corrections to these cross-sections have been computed in |^ for 
ZZ, in I, for W^Z, in |12| for W+W, and in [|l3l Q for W^-f and Z7. These 
computations were all done with the traditional method of evaluating directly the squared 
amplitude through interference (cut) diagrams and evaluating the traces in D = 4 — 2e 
dimensions. As a consequence, the computed cross-sections were summed over all W and Z 
polarization states. 

A more realistic treatment of these processes can be obtained by properly including 
the decay of the vector bosons into massless fermions. In fact, vector bosons are identified 
through these decay products. Thus, the comparison of theory and experiment is much easier, 
since cuts on the kinematics of the decay products can easily be added to the computation. 
As one example of the importance of decay-angle correlations, it has been proposed [^, to 
search for the Higgs boson at the LHC in the intermediate mass range rriH = 155 — 180 GeV 
in the channel H W^W~ i^ui'u. To reduce the continuum W^W~ background, one 
can exploit [|1^] the anti-correlation between W helicities for the signal process, as well as the 
strong correlation between the W helicity and the decay lepton direction; the signal peaks 
when and l~ are nearly coUinear, cos9i+i- ~ 1, while the background is relatively flat in 
this variable. For such a search it is clearly important to understand as well as possible the 
background distributions for cos O^+g^- and other kinematic variables. 

In the narrow-width approximation, vector boson decay is simple to implement at the 
amphtude level. Because the couplings of vector bosons to fermions are spin-dependent 
(especially the purely left-handed W couplings), it is natural to employ the helicity method 
and compute amplitudes for massless external states of definite helicity. The tree-level 
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helicity amplitudes for massive vector-boson pair production with subsequent decay into 



leptons were first computed in [|T^]. The authors of |T3 also showed that the effects of 
decay- angle correlations are significant. 

In [|18], the above two approaches were merged to get a more complete next-to- leading 
order treatment of vector boson pair production. In this work spin correlations were included 
everywhere except for the virtual contribution. Furthermore, the calculations were extended 
to include also non-standard triple-vector-boson couplings. However, only numerical results 
were presented, making it difficult to use these results in future computations. 

The present paper closes this gap by presenting all helicity amplitudes required for next- 
to- leading order (in the strong coupling as) hadronic production of a vector boson pair. In 
particular, we give the one-loop amplitudes with a virtual gluon for the processes 

qq^ W-W+ + il'v'), (1.1) 

qq^ Z Z ^ {a) + (f f), (1.2) 

ud-^ W'Z ^{£u) + {£'£'), (1.3) 

ud W'-f ^ {£y) + 7, (1.4) 

qq^ Z-1 ^ (£l) + 7. (1.5) 

The processes with a Z or a W^'') pair as intermediate state can be obtained from eq. ( |1.3|) 
and eq. ( |1.4| ) by a CP transformation. The decay of the vector bosons into leptons is included 
in the narrow-width approximation. We also present the tree-level amplitudes for the same 
processes with an additional gluon radiated off the quark line. (The gluon may also appear 
in the initial state. The corresponding amplitudes can be obtained by crossing symmetry.) 
Both sets of amplitudes are needed for a complete next-to-leading order computation of the 
cross-sections for vector-boson pair production. 

In the spinor helicity formalism |]19|, the tree-level amplitudes are trivial to obtain and 



the results are very compact. For the one- loop amplitudes, it was not necessary to do a full 
computation. Almost all of the terms could be extracted from one of the helicity amplitudes 
for the process e^e~ — > qqQQ as presented in p^], where q and Q are massless quarks 



of different flavor. In fact, knowledge of the 'primitive amplitude' (see section ^ for the 
subleading-in-color piece of e"'"e~ qqQQ is sufficient to obtain all one-loop amplitudes 
for the processes listed in eqs. ( |1.1|) - (|1.3|) . The amplitude for the process (|1.4|) can also be 



obtained without doing a full computation. It is sufficient to replace the Z in the process ( PTB] ) 
by a virtual photon 7*. Then the desired amplitude can be extracted from the collinear limit 
of the decay products of the virtual photon. Finally, these results allow for the construction 
of the amplitude for the last missing process, eq. (|1.5|). This last result could also be obtained 



from the known subleading-in-color primitive amplitude for e"'"e~ qqg pTl , P2 |. 

Besides their contribution to next-to-leading order pair production rates, the one-loop 
amplitudes also contribute, via their absorptive parts, to kinematic structures which are odd 
under 'naive' time reversal (the reversal of all momentum and spin vectors in a process). 
Such terms are not present at tree level; also, they are washed out if one integrates over all 
the leptonic decay angles. Analogous effects were considered quite some time ago, in the 



production of + 1 jet at hadron colliders p3 



In section ^ we exploit the results of in order to obtain certain primitive amplitudes. 
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called A"" and A^, which serve as building blocks for the construction of all the one-loop 
amplitudes for qq — > W~^W~, ZZ, W^Z —>■ 4 leptons. This construction will be performed 
in section ^. Section ^ is devoted to the processes with a real photon in the final state. 
Finally, in section ^ we compare our results to the literature p, |10|, |12|, |13| and present our 
conclusions. 



2 Primitive amplitudes for VT+iy-, Z Z, W^Z 
2.1 Preliminaries 

It is by now standard to use the helicity method and color ordering of the amplitudes to 



simplify one-loop calculations in QCD (for a review see e.g. |24|). The results for the helicity 



amplitudes will be expressed in terms of spinor inner products, 

{^J)^{k7\kp, [^J]^{kt\kT), (2.1) 

where \kf) is the Weyl spinor for a massless particle with momentum fcj. The spinor inner 
products are antisymmetric and satisfy {ij) [ji] = 2ki ■ kj = Sij. For later use we also define 

m) ^ (Kim;), 

{t\lm\j) ^ {k-\Mm\kj), 

m + m)\3) = {k-Wi+}^^)\kj), (2.2) 
+ rn){n + r)\j) ^ {k^^ +^^){^^ +^,.)\kp, 



and 



Sij — (/uj -|- kj^ , tijm — (^i ~l~ kj -\- A^m) ; 

5i2 = Si2 - S34 - S56, ^34 = S34 - S12 - S56, ^56 = -556 - S12 - S34, (2.3) 

^3 = -512 + -534 + ^56 ~ 2S12S34 — 2S12S56 — 2S34S56. 

For all the processes qq — ^ V1V2 4 leptons, where Vi^2 = or Z, the color ordering of 
the amplitudes is trivial: All diagrams have the same color factor. Thus the full amplitude 
can be obtained by taking the one and only subamplitude (partial amplitude) and multiplying 
it by the overall color factor. 

Still, the diagrams which contribute to gg — > V1V2 — 4 leptons naturally fall into two 
classes, which can be distinguished by their different dependence on the electroweak coupling 
constants, and which therefore are separately gauge invariant. (See Fig. |l].) This decom- 
position is analogous to the decomposition of QCD amplitudes into primitive amplitudes 
p5| . The first class (diagrams (a^) of Fig. |ID has no triple-electroweak- vector-boson vertex. 



There is only one such tree graph, (ao). We call the one-loop box graph (ai) the 'box-parent' 
because the other one-loop graphs in this class, (ai^i) - (ai^3), can be obtained from (ai) 
by sliding the virtual gluon line around, while leaving the electroweak vector bosons fixed. 
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(ao) 



(Vwwv: 



-ivvww: 



-A/wv: 




(bo) 1 (bi) 




Figure 1: (ag) box-parent tree graph; (ai) box-parent one-loop graph; (ai_i) - (ai^s) additional 
one-loop graphs obtained from the box-parent graph; (bo) triangle-parent tree graph; (bi) trian- 
gle-parent one-loop graph. Solid (dashed) lines represent quarks (leptons). In order to distinguish 
the possibly different vector bosons we used zigzag and wavy lines. For the graphs bo and bi we 
show the diagrams for the W^W~ intermediate state. 



(Bubble graphs where a gluon dresses a massless external quark line are zero in dimensional 
regularization.) The tree-level primitive amplitude A^^^*^'"" is defined to be just the contri- 
bution of the graph (ao), omitting all coupling constant prefactors. The one-loop primitive 
amplitude A"' is similarly defined by the sum of graphs (ai) and (ai^i) - (ai^a). The second 
class (diagrams (bj) of Fig. |l|) contains a three-boson vertex. In this class, there are no 
further graphs besides the parent graphs. These primitive amplitudes will be denoted by 
^tree.b j^b ^ Bccausc A^ cousists of a single triangle diagram, it is much simpler than A"'. 

We present amplitudes in the dimensional reduction ||26[ or four- dimensional helicity 
variants of dimensional regularization, which are equivalent at one loop. The conversion to 
other variants is straightforward . 



2.2 Relation to subleading-color e+e qqQQ amplitude 

The one- loop amplitudes for qq V1V2 -^4 leptons can be obtained almost completely 
from the subleading-color amplitude Af {1,2,3,4) for the process e~^e~ qqQQ, as given 
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Figure 2: Parent diagrams for (1,2, 3, 4), i.e. the subleading-in-color part of e+e^ —>■ qqQQ. 
Solid (dotted) lines represent quarks (electrons). Note that the labeling used in this figure and in 
section corresponds to the one used in ref. [^] . 

m eqs. (3.16)-(3.21) of ligS. The 'parent' graphs for that amplitude are shown in Fig. |^. 
Graph (a) is obviously very similar to corresponding box-parent graph (ai) in Fig. ^ The 
quark pair {1,2} can be replaced by a lepton pair, and the virtual gluon that connects the 
two different quark lines can be replaced by an electroweak vector boson. The fact that the 
vector boson can have axial- vector as well as vector couplings, while the gluon has only vector 
couplings, can be accounted for in the helicity formalism simply by dressing the graphs with 
the appropriate left- and right-handed couplings of fermions to electroweak bosons. 

If Ag(l,2,3,4) were truly a primitive amplitude, then one would be able to get qq 
^1^2 — > 4 leptons from it with no further information. However, the Af (1,2,3,4) as given 
in is not truly primitive, for two reasons. 

Firstly, the diagrams (a) and (b) in Fig. ^ are combined together in (1, 2, 3, 4). How- 
ever, it is no problem to remove diagram (b), because its contribution has been identified 
explicitly in |2^ (as the terms in the second bracket in eq. (3.19) for V^^). Both diagrams 
(bi) of Fig. in and (b) of Fig. |^ are QCD vertex corrections, which evaluate to a universal 
factor times the corresponding tree graph. Thus the terms from diagram (b) in Fig. |^, which 
appeared multiplied by Aq'^'^''^^ in [|20|, will reappear here in A^, multiplied instead by ^*ree,6 

Secondly, and more non-trivially, another class of graphs has also been included in 
y4g(l,2,3,4) in addition to those shown in Fig. ^ — those where the quark pair {5,6} 
and the lepton pair {1,2} are exchanged. In the process e~^e~ qqQQ, both classes of 
graphs have equal weight, but in gg — V1V2 ^ 4 leptons they do not (except for the ZZ 
case). The permutation that generates the 'exchange' terms is 

exchange: 1^6, 2 5 . (2-4) 

We need to delete the 'exchange' terms from the (a) terms in Af{l, 2, 3, 4). Since 'exchange' 
takes ti23 to ti24, part of this deletion is easy to implement — we want to omit all terms 
with a cut in the ^124 channel, for example. Also, the second term in the tree amplitude 
^tree.si (3.16) of should clearly be omitted. But this is not enough. There is still 

^In order to simplify the discussion at this point, we use the labehng of in section p.2[ The labehng 
used in the rest of this paper can be obtained from it by the replacement {1, 2, 3, 4, 5, 6} {5, 6, 2, 1, 3, 4}. 
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a potential ambiguity about which terms come from the desired graphs, and which terms 
from the exchange graphs, in the coefficients of ln(si2), ln(s56), ln(s34) and the three-mass 
triangle integral /|™(si2, S34, s^q), as well as the rational-function terms. 

Fortunately, it turns out that all of the ambiguous terms, except for the rational-function 



terms, were given in eqs. (3.20) and (3.21) of ||2^ as the contributions of the desired graphs 



not including the exchange. The exchange that was desired for e^e~ QqQQ was then 
added explicitly, at the end of eq. (3.20). This same statement is not true for a second 
(simpler) form of Af{l, 2, 3, 4), given in eqs. (12.9)-(12.12) of ||22[. In that form, the exchange 
symmetry was used extensively to simplify the coefficients of ln(si2) and /|™, at the cost, 
however, of hopelessly entangling the desired graphs with their exchange. (However, we have 
managed to significantly simplify even the non-exchange form of these coefficients, as one 
may appreciate by comparing eq. ( |2.16| ) below to eq. (3.21) of p][.) 

We have calculated the desired rational-function terms. We find that ^^^'(1, 2, 3, 4) as 
given in eq. (3.20) of EO] is also valid for the non-exchange case, provided only that the 



final '+exchange' instruction is dropped, and that the following rational-function terms are 
added to the existing terms: 

. ^ 1 1 / [16]^ , (25)^ \ _ 1 (2|(1 + 3)|6)(5|(2 + 4)|1) 
2(3|(1 + 2)|4) 1^[12][56] ^ (12)(56)y' 2 ^56 (3|(1 + 2)|4) ' 

2.3 Virtual primitive amplitudes 

After simplifying the coefficients of ln(si2) and /|™, we can write a fairly compact form for 
the required one-loop primitive amplitudes, 

A°(l,2,3,4,5,6), a = a,b. (2.6) 

In the application to the processes gg — V1V2 — 4 leptons, the quark and anti-quark will 
always be drawn from the set {1, 2}, while the leptonic decay products of a vector boson will 
be {3,4} or {5,6}. However, the exact correspondence depends on the fermion helicity; i.e., 
the same function will appear with various permutations of its arguments (1,2,3,4,5,6) in 
section |^. The helicity assignments in eq. (|2.6|) are (1~, 2~^, 3~, 4"*", 5"^, 6~) (when all particles 
are considered outgoing). For the rest of this section the arguments (1,2,3,4,5,6) will be 
implicit. 

First define the fiip symmetry, 

flipi : 1^2, 3^5, 4^6, {ab) ^ [ab] . (2.7) 

The box-parent tree graph is given by 

^,,3,,, ^ ^ (13) [25](6|(2 + 5)|4) _ ^^.g) 

■S34 ^56 tl34 

The triangle-parent tree graph is given by 



A 



tree, 6 



■S12 S34 S56 



(13) [25] (6|(2 + 5)|4) + [24] (16) (3|(1 + 6)|5)1 , (2.9) 
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in a form that agrees with |T] 

At one loop, we perform a decomposition into a universal divergent piece V, and finite 
pieces F", 



cr 



A 



tree, a 



V + iF" 



where a = a,b, the prefactor is 



Cr 



r(l + e)r2(l 



[An] 



2-e 



r(l - 2e) 



:2.10) 



(2.11) 



and we work in dimensional regularization with D = A — 2e. The amplitudes presented here 
are ultraviolet-finite; all poles in e arise from virtual infrared (soft or collinear) singularities. 
The divergent pieces are given by 



3 

Ye 



-Sl2. 
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The finite piece of the triangle-parent contribution vanishes, 

= 0, 

and the corresponding box-parent contribution is given by 

(13)^[25f (2|(5 + 6)14)^(61(2 + 5)|1)^ 

(34) [56] ti34 (1|(5 + 6)|2) [34] (56) ti34 (2|(5 + 6)|1)3 



(2.12) 



:2.13) 



~2mh 

(Si2, ^134! "534, SseJ 



+ 



134 



Li 



-ii34 



2[34](56)(2|(5 + 6)|1) t 



+ 2 



134 



(16) (26) [14]' ti34 Lc 



IM. 



'S34 



[34](56)(2|(5 + 6)|1)2 su 

3 (6|(2 + 5)|4)^ 

4 [34] (56) ti34 (2|(5 + 6)|1) 



(6|1|4)(6|(2 + 5)|4) Lo(^) 
[34](56)(2|(5 + 6)|1) su 

1 (26)[14](6|(2 + 5)|4) 

2 [34] (56) (2|(5 + 6)|1)^ 



In 



"^134) ("-s 12) 



^"■^34 J 



In 



(-^134) (--§12) 
(-^34)' 



J- T 7-3mr, , . ^ I ^ (^234^12 + 2S34S56) f [45] 

+ i J3 l,Si2, S34, Ssej + - 



^34/12 ln( — —) -flipi 

(36)' 



+ 



3 v"iz,^..,"^oD, . 2 (2|(5 + 6)|1)A3 V[34][56] ' (34) (56) 
(36) [45] (ti34 - t234) 1 (6|(2 + 5)|4)2 



:2.14) 



(2|(5 + 6)|1)A3 2 [34] (56) ti34(2|(5 + 6)|1) ' 
where 'flipi' is to be applied only to the terms inside the brackets ([ ]) in which it appears. 
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The coefficient of the logarithm reads 

3 (556 (tl34- t234) (3|(1 + 2)|4)(6|(1 + 2)|5) 
2 



^34/12 



+ 



+ 



+ 



+ 4 



(2|(5 + 6)|l)Ai 
1(3|4|5)[4|(5 + 6)(1 + 2)|5] 
2 
1 



3 (36)[4|(l + 2)(3 + 4)|5] 
^2 (2|(5 + 6)|1)A3 



[56] (2|(5 

i^l34 



2(2|(5 + 6)|1)A3V 



6)|1)A3 

/(34H45f 
[56] 



[14] (26) ti34 ((36) ^i2-2(3|45|6)) 
(56) (2|(5 + 6)|1)2A3 



[34] (36)- 
(56) 



2 (36) [45] 



(3|(1 + 4)|5) (34) [14] (26) \ [45] - 2 [4|36|5] 



(2.15) 



[56] 



(2|(5 + 6)|l)y (2|(5 + 6)|1)A3 



(3|4|5)(6|(1 + 3)|4) + (6|3|4)(3|(2 + 4)|5) 



+ 2 



'12 



(2|(5 + 6)|1)A3 
/[45] (3|(2 + 4)|5) 



(36) (6|(1 + 3)|4) 



(2|(5 + 6)|1)A3V [56] 
and the three-mass triangle coefficient T is given by 
3 S12 S12 (ti34 - t234)(6|(l + 2)|5)(3|(1 + 2)|4) 



(56) 



2 (2|(5 + 6)|l)Ai 

I (3si2 + 2 ti34)(6|(l + 2)|5)(3|(l 
2 



(2.16) 



2)|4) 



+ 



t 



134 



(2|(5 + 6)|1)^A 



(2|(5 + 6)|1)A3 

[14] (26)((3|6|5)556-(3|4|5)534) 



+ 



- [15](23)((6|5|4)556-' 
(36) [45] si2ti34 (34)[56](6|(1 + 2)|4)' 



+ 2 



(2|(5 + 6)|1)A3 
(6|(2 + 5)|4) 



(2|(5 + 6)|1)A3 



(2|(5 + 6)|1)A3 
|5|2)(2|l|4)<556-(( 



+ 2 



3|4)534) 

(16) [24]((6|5|4)556-(6|3|4)534) 



[34] (56) A3 
|2|l)(l|3|4)534+(6|(2 + 5)|4)si2 5i2 



[34] (56) 



2(3|(2 + 6)|5)si2 



[14] (26) (3|(2 + 6)|5_) ^ ^ [15] (23) (6|(2 + 5)|4) 



[14] (26) (6|(2 + 5)|4)5 



12 



+ 



(2|(5 + 6)|1) 
1 



(2|(5 + 6)|1)- 



2(2|(5 + 6)|1) 



'g (6|2|4)(6|l|4) ^ (3|2|5)(3|1|5) 



[34] (56)(2|(5 + 6)|1)^ 
[14] (16) [45] [24] (26) (36) 



+ 



[34] (56) 
(23) [25] (36) 



(34) [56] 
(13) [15] [45] 



+ - 



2(1|(5 + 6)|2) 



(34) 
(16)' [24]' 



[56] 



[34] 
4 (36) [45] 



(56) 



(13)' [25]' 



l [14]'(26)'(ti34^12 + 2 334 356) 

2 [34](56)(2|(5 + 6)|1)^ 



[34] (56) (34) [56] 

For the reader's convenience we recall the definitions of the functions appearing in eq. ( p.l4|) , 

ln(r) , , , _ Lo(r) + 1 



Lo(r) = 



1 — r ' 



Li(r) 



1 — r 
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Ls_i(ri,r2) 



'2mh 



TT 



Li2(l - ri) + Li2(l - rs) + In(ri) ln(r2) - — 

o 



-Uo 1 - 



Li2 1 



(2.17) 



- -In^ 



In 



—m 



In 



-ml 



:2.18] 



where the dilogarithm is 



Li2(x) 



dy 



ln(l - y) 



[2.19) 



The analytic structure of the three-mass-triangle integral /I™ is rather complicated in the 
general case (see Appendix A of ref. [^). However, for the production of a vector boson 
pair V1V2 we have A3 = Si2(si2 - 2Mf - 2M|) + {Mf - M^f > 0, i.e. we are in region (lb) 
of ||20| and the integral is simply given by 

1 



-^3™('5l2! 534, Sse) 



:Re 



2(Li2(-px) + Li2(-P2/)) + In(px) Hpy) (2.20) 



-I- In 



X 



In 



1 + py\ TT 



px 



+ 



where 



X 



S12 S34 

— , y = — , 

■556 -556 



P 



2s 



56 



:2.2V. 



(^56 + 

The rational-function terms in eq. ( p. 141) have been rearranged some, so they no longer 
look exactly like the permutation of the addition of 6p in eq. (|2.5| ) to the rational-function 
terms of eq. (3.20) of [EO . 



2.4 Real (Bremsstrahlung) primitive amplitudes 

A full next-to-leading order calculation of vector-boson pair production requires also the 
tree-level amplitudes with an additional gluon radiated from the quarks. The corresponding 
primitive amplitudes are easily calculated. For the case of a positive-helicity gluon with 
momentum k-j they are given by: 



^tree,a 
^tree,6 



(13) 



(13) [34] [25] (6|(2 + 5)|7) ^ (6|(1 + 3)|4)(1|(2 + 7)|5) 



(17) 534 556^134 

i 

(17) (72) S34 S56tl27 

+ (13)(1|(2 + 7)|4)(6|(3 + 4)|5)- (16) (1|(2 + 7)|5)(3|(5 + 6)|4) 



(2.22) 



^256 (72) 
(36) [45] (1 1 (5 + 6) (2 + 7) 1 1) (2.23) 



The case of a negative-helicity gluon is given simply by applying the operation — flipi to 
the positive-helicity case, where flip^ is defined in eq. ( |2.7] ). These amplitudes do have the 
appropriate limits when the gluon momentum kj becomes soft, or becomes coUinear with a 
quark momentum, ki or k2. 
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3 Dressing with electroweak couplings for 

w+w, zz, w^z 

This section is devoted to the construction of the (squared) amphtudes from the primitive 
amphtudes given in the previous section. We will discuss each process in turn. The formulas 
are given for the tree-level case with only four leptons in the final state. With some trivial 
modifications which will be discussed at the end of this section, the same formulas can be 
used for the one-loop (squared) amplitudes and the tree-level (squared) amplitudes with an 
additional gluon in the final state. 

Concerning the labeling of the external particles: The incoming antiquark always gets 
label 1, while the incoming quark always gets label 2. The external gluon, if present, gets 
label 7. The final-state lepton labelings correspond to the minimal modifications of the WW 
case, 

W- + W+ ^ (4i^4) + {I'^i^',), (3.1) 

namely, 

Z + Z ^ (4I4) + (^^'e), 

W' + z ^ {e^u,) + (44), ^ ■ ' 

and the W^Z case can be obtained by a CP transformation. 

The results will be given in the unphysical configuration where all particles are outgoing. 
The momenta have to satisfy Y^i=\Vi = 0- this configuration the label 1 corresponds to the 
outgoing quark; upon crossing to the physical region it becomes the label for the incoming 
antiquark. Recall that the helicity of the quarks will change sign under the crossing operation. 



3.1 



qq 



WW 



First consider up-quark annihilation into a W^W pair. The leptons (anti-leptons) have to 
be left-handed (right-handed). If the (outgoing) up-quark is left-handed, the tree amplitude 



IS 



tree/ 



A 



where 



c 



S34 



S56 



1 tree,! 



X 



sin 9w ) 

^tree,a^^^ 2, 3, 4, 5, 6) + Cl,„A*'^^^'''(1, 2, 3, 4, 5, 6) 



±2Q sin^ Ow + 



■512 

(1^2Qsin2 Ow) 
^1 



(3.3) 
(3.4) 



si2 - Mi 

6w is the weak mixing angle, and we take (5 = 2/3 and the upper sign in eq. ( |3.4| ) for the up 
quark. The color structure of the amplitude is simply given by 5/^^, where zi, 12 are the color 



labels of the (anti-) quarks. Note that as S12 — > 00, 
cancellation between A^^'^'^'"- and A^^'^^'^: A^^^^'"- 



— > 1, and there is a high-energy 
as S12 — i> C)0. If the up-quark is 
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right-handed, the tree amphtude is 



S34 



■556 



sin^ e 



w , 

1 tree, 6 



S34 -M^ + iTwMw S56 -M^ + iTwMw 



xC^,„A''^<'^'''(2,l,3,4,5,6) 



where 



C 



Ml) 



±2Q sin^ e 



w 



S12 



S12 



(3.5) 



(3.6) 



Since — as as S12 00, the high-energy cancellation is simpler in this case. The 

finite width Vz of the Z boson can safely be neglected in eqs. { WM and ( |3.6|) because we 
have S12 > 4M^. On the other hand, it is of course crucial to keep the iVwMw terms in 
the propagators of the W bosons. 

The corresponding tree amplitudes for down-quark annihilation are 

/ 9 \ 2 



^*-^(rf[,J^;£3,Z/4;f5,^6) 



sin^ e 



S34 



S56 



X 



A'-'''''\l, 2, 6, 5, 4, 3) + Cl4A'''''\1, 2, 6, 5, 4, 3)" 



(3.7) 



and 



^*^'^^(cif,4;£3,^4;^5,^6) 



sin^ e 



X C 



w , 

^tree,6| 



S34 



S56 



S34 -M'+ iVwMw S56 -M'+ iVwM^ 



iw 



2,1,6,5,4,3), 



(3.8) 

where we take Q = —1/3 and the lower sign in eqs. ( |3.4|) and (3^) for the down quark. Note 
that because of Cabibbo-Kobayashi-Maskawa mixing and the large mass of the top quark, 
we have not done the t-channel exchange of the top quark correctly. Fortunately this error is 
proportional to the tiny quantity iVtdl"^ (or for ds annihilation, to VtdV^* times a suppression 
factor for the strange sea quark distribution in the proton). 

We can now construct the differential cross-section in the narrow-width approximation. 
We normalize the squared amplitudes Ai^'^^^ such that the integral over center-of-mass angles 
for the lepton-pair decay products of both vector bosons only has to be multiplied by the two- 
body phase-space factor, in order to obtain the total partonic vector-boson pair production 
cross-section multiplied by the leptonic branching ratios. 



da2{qq ^ V1V2) x Be{Vi) Be{V2) = d^2 J d^^i J ^^^2 W'"' 



The two-body phase-space factor is given by 

/3 



1677 



d COS 9cM , 



where 



P 



'Sl2 - (Ml + M2)2)(S12 - (Ml - M2)2) 



S12 



(3.9) 



(3.10) 



(3.11) 
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and 



dOi sin 9i 

JQ 



2-K 



(3.12) 



Here (^i, <^\) are the angles of one of the decay leptons, measured in that vector boson V^i's 
center-of-mass frame. For equal masses Mi = M2, /3 is the velocity of V\ and V2 in the qq 
center-of-mass frame. 

We obtain for the squared Born amplitude for unpolarized uu WW, 



M^^^'^ui, U2; is, z>4; 4, i^e) = B^W) 



sin^ e 



w . 



w 



(3.13) 



^tree,a^^^ 2, 3, 4, 5, 6) + Cl,uA'''''^\1, 2, 3, 4, 5, 6) + Cr,uA'''''\2, 1, 3, 4, 5, 6) 



while for unpolarized dd WW, we obtain 



M^^^'^du J2; is, 1^4; il i^'e) = Bl{W) 



sin^^ 



w 



Air J 8suNc 



(3.14) 



^^■■-'"(l, 2, 6, 5, 4, 3) + Ci,,A*-^'^(l, 2, 6, 5, 4, 3) + C^,,A*-^'''(2, 1, 6, 5, 4, 3) 



with Nc = 3 for QCD. The corresponding expressions for longitudinally-polarized scattering 
can be obtained simply by dropping the "unwanted" terms in eqs. ( |3.13| ) and ( |3.14| ), and 
adjusting the normalization appropriately. 



3.2 qq^ ZZ ^ Ui'i' 

In the qq —>■ ZZ process, only the symmetric combination appears, 

j^tree,s ^ ^tree,a^^^ 3, 3, 4, 5, 6) + A*'^-.-(l, 2, 6, 5, 4, 3). (3.15) 

(This sum includes the 'exchange' graphs just as in Af; hence the simpler form found for this 
quantity in |22] could be used here, once diagram (b) is removed and the appropriate permu- 
tation is applied.) In the ZZ case, there are more non-vanishing helicity configurations. We 
refrain from writing down all the individual helicity amplitudes in this case, and just present 
the differential cross-section formulas for the various helicity configurations, normalized as 
discussed above: 

M'''"'{qtq^;i^,it;it,ie) = Pzz t;lXe 2, 3, 4, 5, 6)^ , 

M'''%qi:,q^;itri^-Jt,^6) = Pzz<,<e4el^'"''''(l'2,4,3,5,6)p, (3.16) 
M''^%qi:,q^;i^,it;i^,it) = Pzz <,t;iXel^'^''^'^(l,2,3,4,6,5)r, 
M'''''{qtq^;itJ4-J5,^t) = Pzz<X,el^*"'''(l>2,4,3,6,5)p, 



where we defined 



Q N 2 1^-4 
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The formulas for A1*'''^*'(gf , ; . . .) can be obtained by replacing vj^ g by ^ and interchang- 
ing the labels 1 and 2 in eq. (|3.16|) . The unpolarized cross- sect ion, normalized as in eq. ( |3.9|) , 
is given by the sum over all these Ai^^'^'^^s. In the above, the left- and right-handed couplings 
to the Z are 

2 sin^ 9\Y 

(3.18) 



-1 + 2sin^ Ow 

sin 29w ' 
±1 - 2Q sin^ Ow 



sin 26w 
2Q sin^ Ow 



sin 29w ' sin 29w 

where Q is the charge of quark q in units of e, and the two signs in VL,q correspond to up 
(+) and down (— ) type quarks. 



3.3 ud W Z 



Of m 



Finally, for W Z production the down-quark and the lepton I have to be left-handed. The 
lepton which couples to the Z can have either polarization. If it is left-handed, the cross- 
section is proportional to while for a right-handed £' it is proportional to For the 
unpolarized cross-section we have 



W3\ 



3^2 



X < f 



L,e 



Wi,dA*-^'"(l, 2, 3, 4, 5, 6) + i;l,«A*'---'^(1, 2, 6, 5, 4, 3) 

- cot^i^^^^A*'--'''(l,2,3,4,5,6)|' (3.19) 

^^^^^tree,a^^^ 2, 3, 4, 6, 5) + Vl,uA'''''''''{1, 2, 5, 6, 4, 3) 

- cot^w'-^At-^'''(l,2,3,4,6,5) 
The results for the polarized cross-section can easily be extracted from eq. (|3.19| ) 



+ V 



R,e 



3.4 Loop and Bremsstrahlung amplitudes 

The QCD loop corrections to all the amplitudes under consideration are given by a simple 
substitution of loop primitive amplitudes for tree primitive amplitudes. 



^loop 



9 



iV^ - 1 

2 ^ yitree 



a = a,b , 



(3.20) 



where g is the strong coupling, and = 3 for QCD. Thus, the next-to- leading order 
virtual corrections to the above cross-sections are obtained by replacing (^^tree,o-j*^tree,a 



eqs. ( CT) , ( CT) , ( CT) and (^1^) by 2 Re[(A*''<=^'")*A"] and multiplying by ^ - 

The leading high-energy behavior of the loop amplitude cancels in the same way as the 
tree amplitude. First note that none of the terms in F"- contains the factor — - — , which is 



^34^56 ■ 
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present in A^^'^'^'"'. Thus the high-energy hmit is governed completely by the V pieces. This 
ensures that the tree-level cancellation continues to take place at the loop level. 

We remark that any non-Standard Model modifications of the three-boson vertex would 
only afi^ect the second (6) class of diagrams. The virtual QCD corrections can therefore be 
obtained in this case simply by multiplying the modified tree amplitude by the factor V 
given in eq. (|2.12|) . 

The squared amplitudes for the processes with an additional external gluon can be ob- 
tained by a similar change. In this case one has to replace A^'"^''^ and by Aj'"^ and 
^^tree,ay respectively and multiply again by g^^^jf^. Furthermore, su has to be replaced 
by ti27 in the expressions for Cl,Q} and C/j {^}, eqs. (|3^ ) and (|3.6| ), as well as in eq. ( p.l9|) . 
The normalization is the same as long as the two-body phase-space factor d^2 in eq. (|3.10|) 
is replaced by the three-body phase-space factor d^^. Because eqns. ( p. 221) and ( p.23| ) for 
^tree,o crossiug-symmetric, these substitutions are equally valid for processes where the 
gluon is in the final state, qq ViV2g, or in the initial state, gq — > ViV2q and gq ViV2q, 
provided that the factor coming from the spin-color average is changed accordingly; all three 
processes contribute at next-to-leading order. 



4 Real photons in the final state 

The cases with a single real photon, W'j and Z'j, can be handled as well, but the required 
primitive amplitudes are different. The one- loop primitive amplitude needed for the Zj case 
can be extracted from the existing subleading-in-color primitive amplitude for e+e^ qqg, 
as given in appendix IV of This helicity amplitude was first calculated in In the 
W'-y case, just as for the W^W~ and W~Z cases above, the e"'"e~ — >• qqg amplitude has 
to be split further into two 'fully primitive' pieces (which in this case are not related by a 
symmetry). 



4.1 Virtual primitive amplitudes for qq V\'^ 

For qq —>■ V17 there is no 'exchange' symmetry to relate the contributions with the reversed 
ordering of Vi and the photon on the quark line. On the other hand, under an on-shell gauge 
transformation for the photon (i.e., substitution of its polarization vector by its momentum), 
the diagrams where the photon is radiated from the quark line mix with those where it is 
radiated from the Vi = W line. Let the photon be leg 5, the qq pair be {1,2}, and the 
leptonic decay products of Vi be {3,4}. Then we again define two primitive amplitudes, 



A"(l-,2+,3-,4+,5^ 



a 



a, 6, 



(4.1) 



where A'^ (A^) includes all ordered diagrams where the photon is on the leg 2 (leg 1) side of 
the diagram. The negative-helicity photon cases are obtained by another 'fiip' permutation. 



A^(l-,2+,3-,4+,5- 
A^(l-,2+,3-,4+,5- 



fiip2 



A^(l-,2+,3- 



,4+, 5+) 



fiip2U»(l-,2+,3-,4+,5^ 



(4.2) 
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where 

flip2 : 1^2, 3^4, {ah) ^ [ah] . (4.3) 

Therefore we present only the positive-hehcity photon case below. 
The tree amplitudes are given by 



(34) (25) (si2 - S34 



.tree.b _ _ • {'^^f , . 

^ " (34) (15) (.12 - S34) ' ^ ' 

where S34 should be set equal to My^ in accordance with the narrow-width approximation. 

For the one-loop photon amplitudes we use the same decomposition (|2.1CI|) into a divergent 
part V , and finite parts F which are given by 



(13)' , (-S2, , 1 (12)(35)[45][25]Li(zf?t 



(15) (34) (25) -H_,3,'_,3,; + 2 (25) si. 



3 (13)[45]Lo(5ftj 1 [24] [45] 



(25) S12 2 [12] [34] (25) 



(4.6) 



(12)' (35)' 
(34) (15) (25)= 



Ls_i(- 



— S12 — Si5 



■534 "■§34 



1 (23)' [25]' (15) L] 

2 (34) (25) 



^2 
■^15 



(23) [25] ((12) (35) + (13) (25))Lo(5ft 



(34) (25)^ 



Sl5 



1 (12)' (35)' [25] 
'2 (34) (15) (25)' 
1 (23) [25] (35) 



■5l5 



Li(^ 



Lo(ef?f) 



'12 



Sl2 



[14] [25] [45] 



2 [12] (34) (25)' 2 [12] [34] [15] (25) ' 



(4.7) 



The integral functions appearing in the above equations are defined in eqs. ( p.l7| ). 

These amplitudes can be extracted from the collinear limit of the virtual-photon case, 
qq — W^7*, when the momenta of the virtual photon's decay leptons 5 and 6 become parallel. 

The relevant primitive-amplitude combination is A"- H — — 

collinear limit, with fcp = ^5 + fcg, it becomes 



A^ (see eq. (^)). In the 



A'' + 



S12 



-A 



b 5||6 



Sl2 — -534 



^ A^(l-,2+,3-,4+,P+) + L_f A^(l-,2+,3-,4+,p- 



[56] 



(56) 



(4.8) 



Using the 'fiip' relation, we can read off the two independent photon amplitudes. 

The sum of A" and (which is all that is required for the case qq — > Z7) reproduces, 
after relabelings, the existing subleading-in-color primitive amplitude for e^e~ — > qqg, as 



given in appendix IV of [22 
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4.2 Real (Bremsstrahlung) primitive amplitudes 

The tree-level amplitudes with an additional gluon radiated from the quarks, qq —>■ Vi'jg, 
can similarly be obtained from the coUinear limit of the qq Vi'y*g case, which is described 



_^tree,fe^ Let the eluon label be 6. In analogy to the definitions for the 



1 ^ tree, a , 

by ^7 ' + 

virtual case, Aq^'"' {Aq^'^'") includes all ordered diagrams where the photon is on the leg 2 



(leg 1) side of the diagram. The results are 



A 



tree, a/ 
6,7 \ 



l-,2+,3-,4+,5+,6- 



(13)^(21(3 + 4)15) 



AZ''%1-,2 



—I 



^ ,3-,4+,5-,6+) = 
(13) [62] (5|(1 + 3)|4) 

S34 [25] (62) ti34 

-,2+,3-,4+,5+,6+) = 



(34) (25) (16) (62) (ti26 - S34) ' 

(1|(2 + 6)|4)((34)(15) [42] + (35) (16) [62]) 
S34 [25] (16) (62) (ti26 - S34) 



—I 



(13)^(11(3 + 4)15) 



4T''^l-,2^ 



3-,4+,5-,6- 



(15) [24] (3|(1 + 5)|6) 

S34 [15] (16) t234 



(34) (15) (16) (62) (ti26 - S34) 
(l|(2 + 6)|4)((5|(2 + 6)4|3) + (53) ^25) 

S34 [15] (16) (62) (ti26-S34) 



(4.9) 
(4.10) 

(4.11) 
(4.12) 



The cases where the gluon has negative helicity are obtained by applying +flip2 to the above 
amplitudes (which simultaneously reverses the photon helicity). 



4.3 Dressing with electroweak couplings 

As in the case of two massive vector bosons, we present the fully dressed amplitudes and 
differential cross-sections at tree level. The substitutions required to obtain the one-loop 
(squared) amplitudes and tree-level (squared) amplitudes with an additional gluon are ex- 
actly the same as those described in section p.4| . 



4.3.1 ud W~-f 



For the process ud ^ W 7, the lepton (anti-lepton) has to be left-handed (right-handed), 
and the (outgoing) up-quark must also be left-handed. The tree amplitude is 



■S34 



X 



^sin 9w 

g^^tree,a, 



^ a,2,3,4,5^) + QiAf^'''(l,2,3,4,5= 



S34 -Mlr + iTwMw 

■±\ I /O Atrce,h/ 



(4.13) 



where Qi = 2/3 {Q2 = —1/3) is the up (down) quark charge. /^From this amplitude we 
obtain the spin-summed squared amphtude. 



M'"'%u^,d2;h,iy4;i,) 



(4.14) 



ud\ 



sin 9 



w 



J2 \Q2A;{1, 2, 3, 4, 5^) + gi4(l, 2, 3, 4, 5^) 
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The normalization of this cross-section is the same as in the W^W case, eq. ( |3.9[ ), except 
that the integral / (Pfl2 and the factor Bi'{V2) should be omitted. 

4.3.2 qq^ Z-f^ ££j 

In the qq —>■ Zj process, only the symmetric combination appears, 

2, 3, 4, 5±) = ^^=^^^''^(1, 2, 3, 4, 5±) + Af''\l, 2, 3, 4, 5±). (4.15) 

We again refrain from writing down all the individual helicity amplitudes here, and just 
present the differential cross-section formulas for the various helicity configurations, normal- 
ized as discussed above: 

M'''%q^,q^;i^,It;ii) = Pz, Q'viyL^MT''ih2,3,4,5^)\\ 
M'''%qtq^;it,il;it) = Pz, QX, A 2, 4, 3, 5^)^ , (4.16) 

where we defined 

pz^ = B,{Z) eM ^ ) -— ^ ^ . (4.17) 

The formulas for Ai^^^'^{qi, ^2 ; . . .) can be obtained by replacing f £ ^ by f|. ^ and interchang- 
ing the labels 1 and 2 in eq. ( 4.16|) . 



5 Concluding remarks 

We have presented all helicity amplitudes which are needed for a complete computation of 
the next-to-leading order QCD corrections to the production of a Z Z, W^W~, W^Z, W^'~^ 
or Z7 pair at hadron colliders, where the spin correlations are fully taken into account. The 
subsequent decay of each massive vector boson into a lepton pair is included in the narrow- 
width approximation. 

The above cross-sections have been "integrated" over the lepton decay angles, after which 
they reproduce the previously published analytic formulae for the virtual corrections of 
refs. |iy, |l^, |13| to 30 digits accuracy. Such a high accuracy can be achieved because 
the true integration can be replaced by 6 x 6 = 36 numerical evaluations, in which each 
lepton is emitted along three orthogonal axes (both positive and negative directions) in 
the corresponding vector-boson center-of-mass frame. As a check of the bremsstrahlung 



amplitudes we compared in the same way the real contribution in the W W~ case to ref. [0 
and found full agreement. 

Of course the new information provided by the above cross-sections is not the total cross- 
section, but rather the correlations between the lepton decay angles (or lab-frame momenta). 
There are at least two different ways to access this information at next-to-leading order: 

(1) Construct a general purpose Monte Carlo program directly in terms of the lepton mo- 
menta. 

(2) Compute the elements of the density matrix Di-^^i^-j-^j^{s, . . .). These are the amplitude 
interferences for the production of a vector boson Vi with helicity ii or ji = —1, 0, -|-1 along 
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its direction of motion (in the qq —>■ V1V2 frame, for example), and V2 with hehcity ^2 or 
j2- The complete 0{as) density matrix for polarized vector boson production can easily 
be computed using the results presented in this paper. One simply has to carry out the 
integration over the angles of the lepton pair, except that one now weights the numerical 
evaluations with an additional projection operator, exp[z(— +ji)0i + «(— «2 +^2)02], where 
01^2 is the azimuthal angle of a decay lepton with respect to the direction of motion of Vi^2- 
One can then fold the computed production density matrix with the decay density matrices. 
A possible advantage of the second approach is that it is straightforward to include the 
additional QCD corrections that are present for hadronically decaying vector bosons, such 
as those computed in for the W boson. 

For this latter application, to processes such as qq — > W^W" — > q'q" ivi (and neglecting 
interferences with non-resonant processes), the required amplitudes may be obtained from 
those described in this paper by appropriate modifications of the coupling constant factors. 
(At (9(a;s), gluon exchange between the initial and final quark lines gives a vanishing cor- 
rection.) Because quark final states are phenomenologically somewhat less relevant, these 
amplitudes have not been explicitly presented here. Similarly, one can easily extend these 
results to cover the cases where the lepton pair coming from an on-shell Z boson is replaced 
by a Drell-Yan pair of arbitrary invariant mass from an intermediate virtual photon plus Z . 

The amplitudes presented here can easily be implemented in a Monte Carlo program. 
This would extend previous results JlSl in that the spin correlations of the virtual part 
would also be fully included. Furthermore, the simplicity of the amplitudes should aid in 
the next-to-leading-order study of the effects of non-standard triple-gauge-boson couplings. 
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